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The satellite experiment WMAP has produced for the first time a high-coverage, high-resolution 
survey of the microwave sky, releasing publicly available data that are likely to remain unrivalled for 
years to come. Here we compare the WMAP temperature power spectrum, along with an exhaustive 
compilation of previous experiments, to models of dark energy that allow for a tracking epoch at the 
present, deriving updated bounds on the dark energy equation of state and the other cosmological 
parameters. Moreover, we complement the analysis by including a coupling of the dark energy to 
dark matter. The main results are: a) the WMAP data alone constrain the equation of state of 
tracking dark energy to be < — 0.67(— 0.49) to 68%(95%) (confining the analysis to w$ > —1), 

^/-^ I which implies for an inverse power-law potential an exponent a < 0.99(2.08); 6) the dimensionless 

coupling to dark matter is \/3\ < 0.07(0.13). Including the results from the supernovae la further 

f"^ ■ constrains the dark energy equation of state. 
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I. INTRODUCTION 

The WMAP satellite has just released the first full sky maps produced after more than one year of operation 
(Bennett et al. 2003). These data mark an important step in cosmology: they represent in fact the first publicly 
available multi-band high-resolution full-sky survey of the microwave sky. The WMAP survey is likely to remain the 
, highest-quality CMB temperature anisotropy survey until the launch of the Planck satellite not earlier than 2007. 
■ Although the exercise of parameter fitting with each new experiment has produced in the last few years a copious 
literature, the new WMAP data are certainly on a class of its own and deserve to be taken fully into account. Here 
we address the question of what the first year of WMAP data can tell us about the property of dark energy. 

Five years after the first observational hints about the existence of a dominant component of unclustered matter 
L» : . with negative pressure (Riess et al. 1998, Perlmutter et al. 1999), the so-called dark energy or quintessence (Wetterich 
1988; Ratra & Peebles 1988; Frieman et al. 1995; Ferreira & Joyce 1997; Caldwell et al. 1998), there are still very 
. few indications as to its nature. The main reason, perhaps, is that we lack any specific theoretical suggestion on the 
properties of the dark energy, i.e. on its self-interaction potential and on how it interacts with the other cosmological 
components. Several works have tried to constrain dark energy models with pre- WMAP high resolution CMB data, 
leading to bounds on the dark energy equation of state. Unless strong priors on the other cosmological parameters 
are imposed, however, these bounds turned out to be rather weak. In particular, the Hubble constant is strongly 
£2 • degenerated with the dark energy present equation of state, w^, (Huey et al. 1999). For instance, with a flat prior 
h 6 (0.45 — 0.9), Amendola et al. (2002) found that essentially all values < were allowed, while = —0.65 ±0.2 
if h = 0.65 ± 0.05. Similarly, adopting the HST value h = 0.72 ± 0.08 (Freedman et al. 2001) and using updated data 
Melchiorri et al. (2002) found w^, £ (—0.45,-1.81) at 95% c.l. (they included also values < —1). Other works 
that addressed the same issue includes Corasaniti & Copeland (2002), Baccigalupi et al. (2002). With pre- WMAP 
data, only the imposition of a low value of h (e.g. h < 0.7) could allow to distinguish dynamical dark energy from a 
pure cosmological constant (or, more precisely, from a fluid possessing an equation of state close to w<p = — 1 at the 
present). 

In reality, all the constraints for CMB presented in the literature on the dark energy equation of state are necessarily 
somewhat model-dependent, since the fundamental quantity, namely the angular-diameter distance to last scattering, 
in general depends on through two integrals over the cosmic evolution that do not allow to distinguish among 
different parameters along degeneracy curves. Although at low £s the geometric degeneracy could in principle be 
broken by the ISW effect, the cosmic variance and the weakness of the effect make this possibility unrealistic in most 
cases. Here, as in several other works on the topics, we restrict therefore our attention to two classes of models which 
are at the same time general (in the sense of covering most modelizations) and simple. The initial conditions are 
choosen so that the trajectories reach the tracking solutions (notice that off-tracking solutions are indistinguishable 
from a cosmological constant). First, we consider dark energy as a scalar field with inverse power-law potentials 
V ~ <f)~ a (Ratra & Peebles 1988), which recovers the exponential potentials for large values of a ; in this case, the 
equation of state at the present is roughly constant 
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due to the tracking mechanism (Steinhardt et al. 1999). This approximates also the cases in which the dark energy is 
a perfect fluid with constant w^, since the angular size of the acoustic horizon depend on the equation of state only, 
and includes also the cases in which is the f^-weighted average of a slowly- varying function w<p(a) (see e.g. Doran 
et al. 2001). 

The WMAP team analysed their power spectrum including dark energy with a constant equation of state (and 
no coupling). They found for the equation of state < —0.5 with CMB data alone and w^, < —0.78 including 
supernovae (Spergel et al. 2003). Although dark energy with an inverse power-law potential has an equation of state 
which is constant only near the present epoch, we will recover very similar results for this class of models, as expected. 

In the second class we include models of dark energy coupled to dark matter (Amendola 2000; Amendola et al. 
2002). This class of models, in which we include the conformally related Brans-Dicke Lagrangians, has been widely 
studied in the context of dark energy (Uzan 1999, Chiba 1999, Chen & Kamionkowski 1999, Baccigalupi et al. 2000, 
Holden & Wands 2000, Chimento et al. 2000, Billyard & Coley 2000, Bean & Magueijo 2000, Esposito-Farese & 
Polarsky 2001, Sen & Sen 2001, Gasperini et al. 2002, Pietroni 2002) It is important to stress that the behavior of 
coupled dark energy cannot be modeled simply by some choice of the equation of state for <fi because the interaction 
induces on dark matter an effective equation of state different from zero. In this case dark energy mediates a long- 
range scalar interaction that modifies the gravity felt by dark matter particles through a Yukawa-type term (Damour 
et al. 1990). On Newtonian scales, the interaction simply renormalizes Newton's constant for dark matter 

A3 2 

G' = G(l + (2) 

where (3 is the dimensionless coupling, while the baryons remain uncoupled (or very weakly coupled, as required 
by experimental constraints, see e.g. Hagiwara et al. 2002). Then, dark energy effectively violates the equivalence 
principle, but in a way that is locally unobservable. For as concerns us, the main effect of the coupling is to induce an 
exchange of energy between the two dominating dark components so that, after equivalence but before acceleration, 
the universe enters a regime in which not only the equation of state (as in the tracking case) but also the density 
parameters Q mi fl^ are constant. In Amendola (2001) this regime was called 0MDE, to remark the fact that matter is 
not the only dominating component. This epoch ends when the dark energy enters the tracking regime, the potential 
energy takes over and accelerates the expansion. During </>MDE we have O^, = 4/3 2 /3 and an effective equation of 
state for the coupled fluid dark matter/dark energy 

Ptot 4(3 2 

w e = = — 3 

Ptot 9 

Note that this behavior is an example of "early quintessence" (Caldwell et al. 2003), in which the dark energy density 
is not negligible at last scattering. Scope of this paper is to constrain both equations of state, the present one, 
and the past one, w e . This will put limits on the two "fundamental" parameters a, (3. A value a ^ would 
imply that the dark energy is not a pure cosmological constant, while a value (3 ^ would imply a large-scale 
violation of the equivalence principle. In the following we will use interchangeably the parameters w^, w e or a, (3. 
In Fig. 1 we compare the background trajectories for the uncoupled ((3 = 0) and coupled ((3 = 0.1) cases, assuming 
a = l,u>b = 0.02, lo c = 0.1, h = 0.7. Notice the plateau of constant fl^ in the coupled model. 

We start with assuming no coupling, (3 = 0, and varying along with h, u c = n c h 2 ,Ub — fli,h 2 and the slope n s of 
the primordial fluctuations. Then we include (3 and constrain both w<f, and w e . We assume a flat space throughout. 
Notice that our Q c refers to cold dark matter only, so that the total matter content is fl m = fib + £l c . 



II. DATA AND METHOD 



The class of models considered in this paper is the same as in Amendola et al. (2002). We compare the models to the 
combined power spectrum estimated by WMAP (Hinshaw et al. 2003). To derive the likelihood we adopt a version 
of the routine described in Verde et al. (2003), which takes into account all the relevant experimental properties 
(calibration, beam uncertainties, window functions, etc). Since the likelihood routine employs approximations that 
work only for spectra not too far from the data, we run it only for models whose % 2 is less than four times the degrees 
of freedom. We experimented with increasing the \ 2 cut and found no important variations. 

We also compare the data to the whole set of pre- WMAP data. To do this, we use the power spectrum provided by 
Wang, Tegmark and Zaldarriaga (2002), which is a compression of essentially all the available data up to mid-2002 
in 25 £-bins from i = 2 to i = 1700, complete of correlation matrix and window functions (calibration and beam 
uncertainty are included in the correlation matrix). The main entries in this compilation are COBE (Bennett et al. 
1996), Boomerang (Netterfield et al. 2002), DASI (Halverson et al. 2002), Maxima (Lee et al. 2002), CBI (Pearson et 
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Figure 1: Behavior of the energy densities for the coupled (/3 = 0.1, top) and uncoupled ({3 = 0, bottom) model. Notice, in the 
coupled case, the <^>MDE in which Q$ is almost constant and dominated by the kinetic energy of the scalar field. The label T 
denotes the tracking behavior, the label A denotes the final future attractor. 



al. 2002), VSA (Scott et al. 2002). To this we add Archeops (Benoit et al. 2002) with its correlation matrix, window 
functions, beam and calibration uncertainties. For the pre-WMAP data we assumed no reionization (r = 0) because 
this was the best fit before WMAP. 

For the pre-WMAP data we integrated out analytically the calibration and beam uncertainties and the overall 
normalization. Let us denote with N e , a e , a WiC the correlation matrix, the calibration uncertainty and the beam 
uncertainty, respectively, of the e-th experiment, and with I e = {1,1,1,...,} a vector of dimension equal to the 
number of bins in the e-th experiment. Moreover, let C^/ be the theoretical CMB spectrum binned in the £-th bin 
with the experimental window function. Then one gets the remarkably simple likelihood function (see Appendix) 
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where the matrix M e is 



M e = {C'[ d C e , d )- L N e + a\l T e h + o w , e A r A 
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and where Z t = log(C^/) and Zd = logC^ and where the vector A expresses the dependence on £ of the beam 
uncertainty (in the case of a Gaussian beam, A = ...,£f}). This is the likelihood we use for the pre-WMAP 

data. 

Our theoretical model depends on two scalar field parameters, four cosmological parameters and the overall nor- 
malization: 

a,f3,n s ,h,oj b ,ui c ,A. (9) 

As anticipated, to save computing time we found it necessary to restrict the analysis to a fiat space; moreover, we 
fixed the optical depth to r = 0.17, the best fit found by MAP (Spergel et al. 2003). We derived the likelihood also 
for r = 0.1 and found no significative differences for as concerns the two dark energy parameters a, f3. The initial 
conditions for the scalar field are found iteratively for each set of cosmological parameters. The overall normalization 
has been integrated out numerically. Since a, f3 determine uniquely the present and past equations of state, we will 
present the likelihood also as function of w,p 7 w e . We calculate the theoretical Cf jt spectra by a modified parallelized 
CMBFAST (Seljak & Zaldarriaga 1996) code that includes the full set of coupled perturbation equations (see Amendola 
2000 and Amendola & Tocchini-Valentini 2001). The other parameters are set as follows: T cmb = 2.726if, Yhs = 
0.24,^ = 3.04. 

We evaluated the likelihood on two grids of roughly 50, 000 models each (for each normalization): a sparse grid that 
covers a broad volume was used as a preliminary exploration; a second denser grid centered on the peaks of the first 
was then used for the actual calculations. For the first grid we adopted the following top-hat broad priors: (3 G (0, 0.3), 
a G (0.25,20), n s G (0.8,1.2), u b G (0.005,0.04), u c G (0.05,0.3) . For the Hubble constant we adopted the 
top-hat prior h G (0.5,0.9); we also employed the HST result (Freedman et al. 2001) h = 0.72 ±0.08 (Gaussian prior) 
but found only minor differences, since the WMAP results are already very close to the HST distribution. The same 
age constraint (> 10 Gyr) used in most previous analyses is adopted here. 

III. UNCOUPLED TRACKING DARK ENERGY 

We begin by putting (3 — and assuming V((j)) — A(f>~ a . In Fig. 2 we show the likelihood for each parameter, 
marginalizing in turn over the others. The dotted line is for the HST prior h — 0.72 ± 0.08. Notice that, due to the 
degeneracy between and h, the limits on w^, are rather weak, especially if no prior is assumed. The numerical 
results are in Table I. Here and in the following the limits are always at the 68% (95%) c.l. while the errors are at 
68% c.l.. As expected, we find results very close to those in Spergel et al. (2003). The small residual differences are 
probably due to the fact we use a grid instead of a Markov chain and fix r instead of marginalizing over it. 



parameter 


WMAP 


WMAP+HST 


w<t, 


< -0.68(-0.51) 


< -0.73(-0.55) 


a 


< 0.94(1.92) 


< 0.74(1.64) 


h 


69+ - 04 


7n+ - 03 


n s 


1.01±0.022 


1.01±0.022 




0.0247±0.001 


0.0247±0.0008 




19+ 002 


n i 9 +0.015 

u-J^-o.oi 


Table I. Uncoupled dark energy. 



IV. COUPLED DARK ENERGY 

In the flat conformal FRW metric ds 2 = —dt 2 + a 2 dxidx l the scalar coupling modifies the conservation equations 
of dark matter and dark energy as follows: 

4> + 3H<j> + a 2 U^ = C Pc , (10) 
p c + 3H Pc = -Cp c <P, (11) 

where H = a/ a. We assume that the baryons are not directly coupled to the dark energy (otherwise local gravity 
experiment would reveal a fifth force, see Damour et al. 1990; the radiation is automatically uncoupled with this 
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Figure 2: Likelihood functions for uncoupled dark energy models, in arbitrary units. In each panel the other parameters have 
been marginalized. The dotted lines are for the HST prior on the Hubble constant. The horizontal long-dashed lines are the 
confidence levels at 68% and 95%. The vertical long-dashed lines in the panel for wj, mark the upper bounds at 68% and 95% 
confidence levels. 



coupling, see e.g. Amendola 1999). The dimensionless coupling 

* - I5R5- < 12 > 

can be seen as the ratio of the dark energy-dark matter interaction with respect to gravity (Damour & Nordtvedt 
1993, Wetterich 1995). In Amendola (2001) we have shown that the dynamics of the system is insensitive to the sign 
of /3, since the </>MDE and the tracking phases do not depend on it. We will consider only j3 > 0. 

It is important to observe that the CMB bounds on f3 depend on the existence of the 0MDE. This epoch has 
several features that distinguish it from tracking: it is very long (from equivalence to z ~ 10); it depends only on the 
kinetic energy of the dark energy, and therefore is independent of the dark energy potential; it cannot be avoided even 
assuming an extremely small initial dark energy density (contrary to the tracking that may be avoided by selecting 
initial conditions with very low scalar field energy - the so-called "undershooting" trajectories ). Most importantly, as 
shown in Amendola et al. (2002), suffers less from the geometrical degeneracy that plague w,p, since at decoupling 
is non-zero (although there still is a degeneracy (3 — h in the sense that larger h are compensated by larger (3). 
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Figure 3: Likelihood contour plots in the space w^^ trac ki ng ),w e ^ r j > MDE) (or Q,/3) marginalizing over the other parameters at the 
68,95 and 99% c.L. 

In Fig. 3 we show the likelihood for w<f,,w e , marginalizing over the other parameters. We find the following 
constraints at 95% c.L: 

a < 2.08, /3 < 0.13. (13) 

Notice that this implies that an exponential potential (corresponding to a — > oo) is rejected even for (3 ^ 0. In 
place of a and (3 we can use as well as likelihood variables the equation of state during tracking and during 0MDE, 
respectively, using Eqs. Q and (EJ- Then we obtain at the 95% c.L- 

- 1 < W ^tracking) < -0.49, < W e ^ M DE) < 0.0075. (14) 

This shows that the effective equation of state during </>MDE, i.e. between equivalence and tracking, is close to zero 
(as in a pure matter dominated epoch) to less than 1%. The striking difference between the level of the two constraints 
in itm indicates that the CMB is more sensitive to the dark energy coupling than to its potential, as emphasized in 
Amendola et al. (2003). In Fig. 4 we report the likelihood for all the parameters, contrasting the WMAP estimation 
with the pre- WMAP one (the limit on stated in Amendola et al. 2003 was slightly different because here we include 
more pre- WMAP data). As it can be seen, WMAP puts quite stronger limits on w<f, and on (3. In particular, the 
likelihood for now vanishes for w,p > —0.4. 

The other parameters are given in Table II. It appears that the limits on the cosmological parameters n s _u>i,,uj c are 
almost independent of (3, while a non-zero (3 favors higher h. The degeneracy (3 — his reported in Fig. 5. 
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Figure 4: Marginalized likelihood for tracking trajectories. The solid curves are for the WMAP data, the short-dashed curves 
are for the HST prior, and the dotted curves in the panels for w$ and /3 for the pre- WMAP compilation . The horizontal 
long-dashed lines are the confidence levels at 68% and 95%. The vertical long-dashed lines in the panel for mark the upper 
bounds at 68% and 95% confidence levels. 



parameter 


WMAP 


WMAP+HST 


pre- WMAP 


w^, 


-0.67(-0.49) 


< -0.69(-0.52) 


< -0.50(-0.25) 


a 


< 0.99(2.08) 


< 0.90(1.84) 


< 2.0(6.0) 


W e 


< 0.0025(0.0075) 


< 0.0023(0.0075) 


< 0.0075(0.016) 


P 


< 0.075(0.13) 


< 0.072(0.13) 


< 0.13(0.19) 


h 


0.73± 0.05 


0.73± 0.04 


>0.62(0.55) 


n s 


1.019± 0.025 


1.018±0.025 


0.97±0.03 




0.0247± 0.0008 


0.0250±0.0008 


0.021±0.003 




0.123± 0.016 


0.120± 0.016 


0.12±0.04 


Table II. Coupled dark energy. 



In Fig. 6 we show the contour plot of the likelihood L(Q m , w<f,), where Cl m = fl c + fib, marginalizing over the other 
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Figure 5: Likelihood for f3, h. This shows the residual degeneracy between the two parameters due to the geometric degeneracy 
in the angular diameter distance to last scattering. 

parameters. We also add the confidence region from the Hubble diagram of SNIa (we used the fit C of Perlmutter et 
al. 1999, plus the supernova SN1997ff at z w 1.75 , Benitez et al. 2002). The product of the two likelihood functions 
is shown in the same figure. It turns out that fl^ = 0.67 ± 0.05 and < —0.76 (95%c.l.) (see Fig. 7). The limit on 
w<p is very close to that obtained in Spergel et al. (2003): this shows that this bound is almost independent of the 
coupling (3. 



CMB observations are a powerful probe to constrain the properties of dark energy. In particular, since a dark 
matter-dark energy interaction would obviously escape any local gravity experiment, cosmological observations like 
the CMB are the only way to observe such a phenomenon. Since observations require the baryons to be decoupled 
from dark energy (or coupled much more weakly than dark matter) , the search for a non-zero coupling (3 is also a test 
of the equivalence principle. We found that current CMB data are capable to put an interesting upper bound to the 
present dark energy equation of state (all limits to 95% c.l.) 



which becomes < —0.76 taking into account SNIa data (let us remember that in our tracking models is confined 
to be > — 1) and to the total dark energy density, 



V. 



CONCLUSIONS 



W0 < -0.49, 



(15) 



fi = 0.67 ±0.05. 



(16) 



For the dark matter - dark energy coupling we obtain: 



\f3\ < 0.13, 



(17) 



regardless of the potential (within the tracking class we considered) , corresponding to a past equation of state < 
w e (<j)MDE) < 0.0075. This implies that the scalar gravity is at least 1//3 2 w 60 times weaker than ordinary tensor 
gravity. As shown in Amendola (1999), the limit on (3 can be restated as a limit on the constant £ of the non-minimally 
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Figure 6: Contour plot of the likelihood function L{Q, m ,w,f,). The dotted lines are for WMAP only, the thik gray lines for the 
supernovae la, and the gray regions are for the combined WMAP+SNIa (all contours at the 68% and 95% confidence levels). 



coupled gravity, £ < 0.01. We have shown in Amendola et al. (2003) that an experiment like the Planck mission can 
lower the upper bound to (3 to 0.05. a limit comparable to the constraint that local gravity experiments impose on 
the scalar gravity coupling to baryons, 0l aryons < 10~ 3 (see e.g. Hagiwara et al. 2002). 
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Appendix 

Suppose we have the binned observational spectra C' td = (l+b e )Ci t d and the binned theoretical spectra C' t t = btCij 
from a given experiment labelled e, where b e is a calibration factor which is assumed to be distributed as a Gaussian 
with zero mean and variance a e and b± is the overall theoretical normalization, which is uniformly distributed in 
(0,oo). Let N e be the correlation matrix for the e-th experiment and 

K = (cjcy- 1 ^ (is) 

the approximated correlation function for the lognormal variables Z = log C (here and in the following we suppress 
for clarity the indexes W of the correlation matrices and the index I from the vectors Z,Ct,Cd) ■ If the assumed 
likelihood is log-normal then we can integrate out b e as follows 

L e = f \{db e e-^) 2 x (19) 

J e 

CX P I - \ E W - b -- Z ^ N ^ W -be-Z d ]\= f db t [] L e (20) 
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Figure 7: Likelihood for fi m and w$ including the constraints from supernovae la (marginalizing over all the other parameters, 
including (3). 



where Z[ = log b t C t and b e « log(6 e + 1), Z& — logC^ (for small b e ) so that 

1 



exp 



(Z' t - z d ) w- 1 {Z' t - Z d ) 



(21) 



where 



iy e = (CjC d )- 1 iV e + a e 2 / T / 



where / = (1, 1, 1, n e ), n e being the number of datapoints in the e-th experiment. Note in fact that (Woodbury 
formula, see e.g. Melchiorri et al. 2002) 

(N + XX T ) = N- 1 - [7Y- 1 Z(1 + X T N- 1 X)- 1 X T N~ 1 ] 

which gives the above if X — a e I. 

Suppose now the experimental spectra are distorted by a further ^-dependent factor (l+w e Ai) due to the uncertainty 
on the beam angular size, where w is a random variable distributed as a Gaussian with variance a w ^ e and Ag is a 
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constant vector determined by the experiment. Then each likelihood L e should be further integrated as follows 



L' 



dw P e 



exp 



(w e A + Z' t - Z d ) W~ l (w e A + Z' t - Z d ) 



where, on applying again the Woodbury formula, 



L' e = exp 



--{Z' t -Z d )M-\Z[-Z d ) 



where 



M e = (CFC)- l N e + a 2 J T I + a 2 w e A T A 



(22) 

(23) 
(24) 



(25) 



(26) 



In the case of a Gaussian beam of angular size 9q the spectrum is multiplied by a factor e l e o . It can be seen then that 
a small misestimate 9 (assumed to be distributed as a Gaussian variable with variance ag ) of the beam size induces 
a correction on the spectrum at the £ — th bin by a factor [1 + £ 2 (9\ — 2 )] = 1 + w£ 2 , where the random variable 
w « 28o(do — 0) is also Gaussian and has uncertainty a w — 29 ag. The vector At can be approximated then as £ 2 . 
Finally, we integrate over b t = \ogb t . Then we have Z t = b t + Z[ and 



L = J db t ~[[L' e (bt) = J db t exp -Ij2(bt + Z t - Z d ) M' 1 (b t + Z t - Z d ) 



exp 



° 1 1 a 



(27) 

(28) 
(29) 



where 

a = 

e 

j3 = J2l T M-\Z t - Z d ) 

e 

e 

which is the formula we use in this paper for the pre-WMAP data. 
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